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. Abstract 

• Algorithms are described that help with obtaining a classification of the semisim- 

, pie subalgebras of a given semisimple Lie algebra, up to linear equivalence. The 

algorithms have been used to obtain classifications of the semisimple subalgebras of 

the simple Lie algebras of ranks < 8. These have been made available as a database 

inside the SLA package of GAP4. The subalgebras in this database are explicitly 

' given, as well as the inclusion relations among them. 

(N ■ 

1 Introduction 

O I It is an extensively studied problem to classify the semisimple subalgebras of a complex 
semisimple Lie algebra g, up to an equivalence relation. The most natural equivalence 
relation for this is the one of conjugacy by the inner automorphism group G. Two sub- 
algebras of g are simply called equivalent if they are conjugate under G. In [5], Dynkin 
also considered the notion of linear equivalence: two subalgebras gi, g2 C g are said to be 
linearly equivalent if for every representation p : q ^ qK^) subalgebras p(gi), p{Q2) of 
g[(y) are conjugate under GL{V). 

A subalgebra of g is called regular if it is normalised by a Cartan subalgebra of g. 
Semisimple subalgebras of this kind correspond to root subsystems of the root system of 
g. An S-subalgebra is a subalgebra which is not contained in a regular subalgebra. 

In |1] Dynkin classified the maximal semisimple 5'-subalgebras of the Lie algebras of 
classical type, upto equivalence. More precisely, P] contains a description, or a procedure, 
by which it is possible for a given Lie algebra of classical type to find its maximal semisimple 
S'-subalgebras. 

Dynkin treated the Lie algebras of exceptional type in [5]. The main results of this 
paper are 

• an algorithm to classify the regular subalgebras of a semisimple Lie algebra, up to 
equivalence. 
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• a classification of tlie semisimple S'-subalgebras, up to equivalence, of the Lie algebras 
of exceptional type, 

• a classification of the simple subalgebras, up to linear equivalence, of the Lie algebras 
of exceptional type. 

Lorente and Gruber ([I2]) applied Dynkin's methods to obtain explicit lists of semisim- 
ple subalgebras of the simple Lie algebras of classical type. More in particular, they 
obtained lists of the regular subalgebras and of the S'-subalgebras of the Lie algebras of 
classical type of ranks < 6. 

Recently Minchenko ([H]) has revisited Dynkin's classification of the simple subalgebras 
of the simple Lie algebras of exceptional type. He corrected several small mistakes (most 
notably he found two extra simple subalgebras in the Lie algebra of type Eg). Secondly 
he found the classification of the simple subalgebras up to equivalence. Thirdly, he has 
computed a lot of additional data (such as the normalisers of the subalgebras in G). 

One motivation for studying semisimple subalgebras of semisimple Lie algebras comes 
from theoretical physics. In models like the vibron model and the interacting boson model 
(cf. |10]) chains of subalgebras are used. For applications of this kind the subalgebras need 
to be explicitly given, i.e., for each equivalence class a representative needs to be given by 
a basis. Furthermore, methods are needed to obtain the inclusion relations between the 
subalgebras (more precisely: to decide whether two given classes have representatives such 
that one is contained in the other). The classifications present in the literature do not 
appear to immediately give this. For example, in [5], only the S'-subalgebras are explicitly 
constructed. And only the simple subalgebras are listed, and not the semisimple ones (with 
the exception of the S'-subalgebras). Finally no inclusion relations are given (again with 
the exception of the S-subalgebras). 

The aim of this paper is to describe algorithms, and report on the results obtained 
with their implementation, that help with obtaining a classification of the semisimple 
subalgebras of a given semisimple Lie algebra, up to linear equivalence. Furthermore, the 
subalgebras are explicitly constructed, as well as the inclusion relations among them. Here 
we say that the algorithms "help" to obtain a classification as one step in the algorithms 
(the construction of the subalgebras) is not entirely algorithmic - occasionally some human 
intervention is needed for that. 

Equivalence implies linear equivalence, but the converse is not always true. However, 
if Q is of type An, Bn, Cn, F4, G2 then the two concepts coincide (cf. [H], Theorem 3). 
In the remaining types there are some exceptions and they are explicitly described (|5], 
|14]). Hence it is straightforward to obtain the classification of the semisimple subalgebras 
up to equivalence from the list of semisimple subalgebras up to linear equivalence. One 
of the main advantages of linear equivalence as opposed to equivalence is that we have a 
method for deciding it (see Section [3]). For these reasons in this paper we focus exclusively 
on linear equivalence. 

By considering embeddings of Lie algebras in g, rather than subalgebras of q we get a 
slightly different perspective on the problem. Also for embeddings we have the notions of 
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equivalence and linear equivalence. Let g be a semisimple Lie algebra, and (pi,y^2 '■ Q "-^ Q 
injective liomomorpliisms. They are said to be equivalent if there is a a G G with tpi = aLp2- 
They are said to be linearly equivalent if for each representation p : g — 0l(^) the induced 
representations pipi, pip2 of q are equivalent. Let g be a semisimple Lie algebra, and g' C g 
a subalgebra isomorphic to g. There can be several non-equivalent embeddings g g'. 
This is only possible if g has outer automorphisms. From a classification of subalgebras up 
to linear equivalence it is straightforward to get all embeddings up to linear equivalence. 
For this reason we concentrate on constructing subalgebras, rather than embeddings. 

One approach to the problem is to start from the existing classifications in the lit- 
erature. One could take the maximal S'-subalgebras constructed by Dynkin, along with 
the regular subalgebras, and by successively constructing their subalgebras get the entire 
list of subalgebras. This, however, would not confirm, or correct, the existing classifcation. 
Moreover, if the list of maximal subalgebras has an error, then this will lead to many errors 
in the resulting classification. (And it appears that this can, for example, easily happen 
in type D2ni see below.) For these reasons the approach taken here aims at obtaining the 
classification from scratch. This has the added advantage that the classifications in the 
literature and the new ones can validate each other. In particular, if both are the same 
then this constitutes a good argument for their correctness. 

The main idea used here to classify subalgebras is to start with the ones of smallest 
rank. The subalgebras of rank 1 are well-known from the classification of the nilpotent 
G-orbits in g. Secondly we construct the subalgebras of higher rank as a kind of extension 
of algebras of lower rank. This way we "climb our way up". So, in a sense, it is the reverse 
approach to starting with the maximal subalgebras. 

The algorithms described in this paper have been implemented in the language of the 
computer algebra system GAP4 ([6j), using the package SLA The main result that 

has been obtained using this implementation is a database of all semisimple subalgebras 
of the simple Lie algebras of ranks < 8. This database is also contained in the package 
SLA. It also contains all inclusion relations between the linear equivalence classes. It is 
complemented by a function for computing the semisimple subalgebras of a semisimple, 
non-simple. Lie algebra. In Table [T] we show some statistics relative to the simple Lie 
algebras of ranks 7, 8. The table contains the number of (linear equivalence classes of) 
subalgebras, and the number of their isomorphism types. 

There is also the question of the field of definition. The simple Lie algebras are given 
by a multiplication table relative to a Chevalley basis. The subalgebras are given by a 
basis. However, not all linear equivalence classes of subalgebras have a representative with 
a basis with coefficients in Q (with respect to the given Chevalley basis). Our results show 
that for all semisimple subalgebras of the simple Lie algebras of ranks < 8 there exists a 
field extension F of degree < 2 of Q, such that the subalgebra can be given by a basis 
with coefficients in F. The last column of Table [1] gives a field extension F of Q such that 
all semisimple subalgebras can be given by a basis with coefficients in F. Here we remark 
that it is by no means clear that these are the smallest possible fields (except, of course, 
when the field is Q). We have made an effort to keep the fields small; but the problem of 
finding the absolute smallest field is a difficult one which we do not solve here. 
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Table 1: Semisimple subalgebras of the simple Lie algebras of rank 7 and 8. The second 
column displays the number of linear equivalence classes of subalgebras. The third column 
has the number of different isomorphism types of subalgebras. The last column has a field 
over which all subalgebras are simultaneously defined. 



Next there is the question of the validation of the results: how can we be certain that our 
classifications are correct? Although in this paper we prove the correctness of the method 
that we use, there is still ample possibility to make mistakes while using it. However we 
do have some circumstantial evidence for the correctness of our lists. Firstly, the method 
does not deal with regular subalgebras differently than with other subalgebras. But at the 
end we find the same regular subalgebras as with Dynkin's algorithm from [5j. Secondly, 
the S'-subalgebras that we find in the exceptional types coincide with the ones found by 
Dynkin. In the classical types for ranks < 6 we find the same S'-subalgebras as Lorente and 
Gruber ([12]) (except in D4, De, see below). Thirdly, also the lists of simple subalgebras 
agree with those found by Dynkin (and in the case of £^8 corrected by Minchenko). 

One result of our calculations is that in type -D2n, for n = 2,3,4, there appear maximal 
semisimple subalgebras which are isomorphic, but not linearly equivalent. In type there 
are three (linear equivalence classes of) maximal subalgebras of types AxB^ and Bj,. In 
type there are two maximal subalgebras of each type AiCa and A^. And in there 
are two maximal subalgebras of each type B2B2, B4, A1C4 and A^. This appears not to 
have been known in the literature, for example [13] lists one algebra for each of the above 
types. In all cases the algebras are conjugate under outer automorphisms. It would be 
interesting to formulate and prove a general statement about the maximal subalgebras of 
the Lie algebra of type D2n- However, this would be beyond the scope of this paper. We 
intend to come back to it in a subsequent paper. 

This paper is structured as follows. The next section describes a number of concepts 
and results from the literature that we need. This allows us at the end of the same section 
to give a more or less detailed description of the method we use. The subsequent sections 
then describe every step in detail. 

Acknowledgement: It is my pleasure to thank Luigi Scorzato and Lorenzo Fortunato 
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for having suggested the topic of the paper to me, and for many fruitful discussions about 
it. 



2 Preliminaries 

Throughout Q will be a semisimple Lie algebra over C, with a fixed Cartan subalgebra f). 
The inner automorphism group of g will be denoted G. 

2.1 The action of the Weyl group 

The Killing form on g will be denoted k; it is defined by k{x, y) = Tr(ada; • ady). The form 
K, is nondegenerate on q and on f). Hence we can define a bijection f)* — > f), i->- /i, by 
K{fi, h) — ii{h), for /i e f). Then (//, A) = A) defines a non-degenerate bilinear form on 
f)*. Also for 7^ we set 

For a e f)* we set 0a = {x G | [h, x] = a{h)x for all h G i)}. We let $ be the set of all 
nonzero a e f)* with 0a 7^ 0. Let f)^ be the real vector space spanned by Then ( , ) is 
an inner product in 1^^, and $ is a (reduced) root system in f)^. 
For a,P ei)^ we set 

For a G $ we define the refiection Sa ^ f)^ — ^ by Sa(/u) = ^ — {^,a'^)a. The group 
generated by all Sa for a G $ is called the Weyl group, and denoted W. 

Let l)]R be the real vector space spanned by all a^, for a G For a G $ we define the 
linear map Sq : f)]R — >■ f)]R by Sa{h) — h — a{h)ha- Then SaiP^) — — (a, /3^)q;^. A small 
calculation shows that the following diagram commutes 

Or ^ Or 



Or ^ Or- 

So, more generally, for w G and G f)^ we have wfi — wjl. Also leaves the Killing 
form on f)^ and on f)^ invariant. 
Let < be an order on f)^ with 

• u < V implies u + w < v + w ior all w & f)^, 

• u > implies Xu > for all positive A G M, and Xu < for all negative A G M. 

We call such a < a root-order. A root-order defines a partition $ = U $~ of $ into 
positive and negative roots, and a set A of simple roots. Conversely, if A is a set of simple 
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roots, then we can define an ordering as follows: express u, v as linear combinations of the 
elements of A, and set m < f if the first nonzero coordinate of — m is positive. This order 
will then yield A as set of simple roots. 

Let C* C f)R be set of all /i with (/x, a"^) > for all a > 0. Then every W^-orbit in {)^ 
has a unique point in C*. It is called the fundamental Weyl chamber of f)^. Also we let 
C C l)M.he the set of all h with a{h) > for all a > 0. Again, every ly-orbit in has a 
unique point in C. It is called the fundamental Weyl chamber in l)^. 

2.2 Nilpotent orbits 

Let e G be nilpotent; then the orbit G ■ e is called a nilpotent orbit. Here we recall some 
facts on the classification of nilpotent orbits from [1], [2]. 

Let e G be nilpotent, then by the Jacobson-Morozov lemma there are /i, / G g with 
[h, e] = 2e, [h, f] = —2/, [e, /] = h. We say that {h, e, /) is an s(2-triple. Note that G 
acts on s[2-triples by a ■ {h,e, f) = {a ■ h, a ■ e, a ■ f). Let e, e' G g be nilpotent, lying in 
s[2-triples {h,e,f) and {h',e',f'). Then the following are equivalent: 

• e,e' lie in the same G-orbit, 

• {h, e, f) and {h', e', /') lie in the same G-orbit, 

• h, h' lie in the same G-orbit. 

Let {h, e, /) be an 5[2-triple in q. Then h lies in a Cartan subalgebra of q. As all Cartan 
subalgebras of q are G-conjugate, after possibly replacing the triple by a G-conjugate, we 
may assume that /i G P). Then /i G f)M (indeed: a{h) G Z for all a G $). Two elements 
of l)iR are G-conjugate if and only if they are VT-conjugate (cf. [2] Theorem 2.2.4). Hence, 
after a further conjugation we may assume /i G G. In fact, this h determines the nilpotent 
orbit uniquely; it is called the characteristic of the orbit. 

We call an /i G () admissible if it lies in an s[2-triple {h, e, /). Let ei, . . . , be rep- 
resentatives of the nilpotent G-orbits in g, lying in sl2-triples {hi,ei,fi), with hi G G. 
Then 

t 

n = [jw-hi 

i=l 

is the set of all admissible elements in f). 

We will often have the need to run through a 1^-orbit W ■ h^. For this Snow ([16], 
see also [7]) has devised an efficient algorithm, which makes it possible to run through the 
orbit and inspect each element without storing all of the orbit. This feature will be very 
important for us. 

2.3 The Dynkin index 

Assume that g is simple. It is well-known that upto multilication by nonzero scalars, there 
exists a unique nondegenerate symmetric G-invariant bilinear form on q. The Killing form 
is such a form. 
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Let g C 5 be a simple subalgebra. Let G be the group of inner automorphisms of g. 
Then G G G. Hence the Kilhng form k of g induces a G-invariant bihnear form on g. Let 
R denote the Kilhng form of g. So k{x, y) = 7]k(x, y) for ah x,y E g, where t] is a nonzero 
scalar. 

If we normalise k so that ^(q!'^,^'^) = 2 for the short roots a, and do the same for 
K, then r] is called the Dynkin index of g in g. It is the same for all G-conjugates of g. 
However, it can also happen that nonconjugate subalgebras have the same Dynkin index. 

Lemma 1 Let g C g be a semisimple subalgebra that is the direct sum of simple ideals, 
g = gi ® ■ ■ ■ ® Qm- Then K{gi,gj) =Oforiy^j. 

Proof. Let x,y G g^ and z G g^. Then k{[x, y], z) = k{x, [y, z]) = 0. So since gj = [g^, g^] 
the result follows. □ 



2.4 Canonical generators 

Let g be a semisimple Lie algebra. Then g has a canonical set of generators ([llj, Chapter 
IV). That is a a set of elements xi, . . . ,Xr, yi, . . . ,yr, hi, . . . ,hr such that 

[hi, hj] = 

[xi,yj] = 6ijhi (1) 
[hj,Xi] = C{i,j)xi 
[hj^TJi] = -C{i,j)yi. 

Here G is the Cartan matrix of the root system of g. We call the sequence {hi, . . . , hj.) the 
/i-part of the canonical generating set. We note that hi = , where {ai, . . . , a^} is a set 
of simple roots of the root system of g. 

Suppose now that /ij G f). Then the /ij G f) are admissible; hence lie in the set "H of 
Section 12.21 In the sequel we will say that the /i-part of a canonical generating set lies 
in i) to mean that all of its elements do. The next theorem is essentially the same as 
Theorem 4.1. 

Theorem 2 Let G be the Cartan matrix of a root system. Let Xi, yi, hi be elements of a 
finite dimensional Lie algebra satisfying the relations (QP- Then the subalgebra generated 
by these elements is semisimple, and its root system has Cartan matrix equal to G . 

Proof. For i ^ j consider the element 

y,, = (ad^.)-^(^'^)+^fe). 

A short calculation (cf. [7], Lemma 7.n.3) shows that [xi, yij] = and [hi, yij] = {G{j, i) — 
'2)yij. But G{j,i) — 2 < 0. It follows that yij generates a finite-dimensional irreducible 
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st2-niodule of negative highest weight. This is impossible, hence yij = 0. Similarly we 
have 

(adXi)-^(-''')+i(x,) = 0. 

Hence the Xi, fji, hi satisfy the Serre relations (see [12], Chapter VI, §4). This implies that 
the algebra they generate is a quotient of the semisimple Lie algebra u corresponding to 
the Cartan matrix C by an ideal. This ideal is the sum of some of the simple ideals of u. 
But since the Xi, jji, hi are nonzero, this ideal has to be zero. □ 

2.5 Solving polynomial equations 

In order to construct the subalgebras that we are after, on some occasions we need to solve 
polynomial equations in several variables (see Section H]). For this no general algorithm 
exists, so we have to do it by hand. However, a computational tool that makes this a lot 
easier is provided by Grobner bases. 

Let /i, . . . ,/s € k[xi, . . . ,Xn], where is a field, generate the ideal /. Then solving 
/i = • ■ ■ = = is the same as solving g = for all g G Q, where Q is any other 
generating set of /. A Grobner basis is, on many occasions, a particularly convenient 
generating set for this purpose. Especially if the Grobner basis Q is computed relative to a 
lexicographical ordering, then Q has a triangular structure, which often makes solving the 
polynomial equations easier. Also, if there are no solutions over the algebraic closure of k, 
then the reduced Grobner basis is {1}. So this situation is immediately detected. Here we 
do not go into the details, but refer to [3] for an in-depth discussion of Grobner bases and 
polynomial system solving. 

2.6 Outline of the method 

Here we summarise the method we use to classify semisimple subalgebras of g. 

Let C be the r x r Cartan matrix of the root system of a semisimple Lie algebra. The 
objective is to classify the semisimple subalgebras of q having a root system with Cartan 
matrix C, up to linear equivalence. We assume that the semisimple subalgebras of g of 
smaller rank have been classified. We note that the classification for rank 1 is known from 
the classification of the nilpotent orbits in g. 

Let g C be a subalgebra with Cartan matrix C, and canonical set of generators hi, 
Xi, Vi, i < i < r, satisfying ([1]). Then the hi lie in a Cartan subalgebra of g. So since 
all Cartan subalgebras of g are conjugate under G, we get that g is equivalent, and hence 
linearly equivalent, to a subalgebra with a canonical generating set with the /i-part lying 
in f). So we may assume that /li G f), and hence hi eH. 

In Section [5] we describe methods to assemble a set H of r-tuples {hi, ... , hr) G 'H'' 
such that all classes of linearly equivalent subalgebras with Cartan matrix C have a repre- 
sentative that has a canonical generating set with /i-part in H. Here one of the objectives 
is to keep this set "small" . 
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Let {hi, . . . , hr) E H. Section H] contains methods that construct Xj, i/i in g satisfying 
(II]), or decide that no such elements exist. In the former case we have found a semisimple 
subalgebra of g with Cartan matrix C by Theorem [2l In the latter case the hi do not form 
the /i-part of a canonical generating set of a subalgebra with Cartan matrix C. 

In Section |3] we describe a method for deciding whether two semisimple subalgebras 
are linearly equivalent. So we can get rid of any linearly equivalent pairs of subalgebras 
constructed in the previous step. In fact, linear equivalence depends only on the /i-parts of 
the canonical generating sets; so we can construct the set H so that no linearly equivalent 
subalgebras arise. This is important as constructing the subalgebras is one of the most 
difficult steps. 

We used these methods for classifying the semisimple subalgebras of the simple Lie 
algebras of ranks up to 8. For classifying the semisimple subalgebras of the semisimple, 
but not simple. Lie algebras we have a separate method, described in Section [61 Finally 
the last section has the algorithm that we use for deciding inclusion. 

3 Deciding linear equivalence 

The purpose of this section is to describe an algorithm for deciding whether two semisimple 
subalgebras of g are linearly equivalent. For this we assume that they are given by canonical 
sets of generators, with the /i-parts lying in f). First we prove a theorem that in essence is 
due to Dynkin ([5], Theorem 1.5). Here we show how Dynkin's argument can be adapted 
to prove the statement that we need (Corollary [5]). For this the language of embeddings is 
more appropriate. 

Let : g — )• g be an embedding of the semisimple Lie algebra g into g. Let [) he a. fixed 
Cartan subalgebra of g and assume C f). Let t)^, t)^ be the M-span of the roots of g 
and g respectively. We define a map : {)^ — )■ f)j^ by Lp*{n){K) = fi{ip{h)), where h G ^j^. 

Let pig—)- gt(f/) be a representation and let /i be a weight of p, with weight vector 
V. Then for /i G f)^ we get p{ip{h))v = p{if{K))v = Lp*{fi){h)v. It follows that f*{p) is a 
weight of the representation p^p of g. In particular, it lies in ^j^. Since the weights span 
the spaces f)^, f)^, it follows that 'p*{i)^) = f)^. 

Lemma 3 Let the notation be as above. Fix a root-oder -< o/ flj^- Fix also a division 
$ = $+ u $^ of the roots ^ of g into positive and negative roots, corresponding to a root- 
order of l)^. Let A C $ denote the corresponding set of simple roots. Then there exists a 
root-order < of and a a E Ng{^) with the following properties: 

• The set of positive roots with respect to < is also $+; 

• for tf) = (Tip we have that ^*{p) -< implies /i < A, 

• if tp*{pi) 7^ ip*{\) then p, < X implies tp*{p) -< '0*(A). 

Proof. Let <' be any root-order on f)j^. Define the root-order <" on f]J by p <" X if 
</j*(/i) -< (/?*(A), or if those two are equal, p <' X. Let A" be the corresponding set of simple 
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roots. Then there is w & W with wA" = A. Let a G Ndi)) be such that the restriction of 
cr to [) is w. Set = <^'^, and define the root-order < by: /i < A if <" w^^X. 

Let a E A, and write a = wP for some /3 G A". Then w~^a = j3 >" 0. Hence a > 0; 
and therefore the set of positive roots with respect to < is 

Next, using {w~^fj,){h) = fx{wh) for /i G f) we get ip*{fJ.) = Lp*{w^^fi). Hence ip*{fi) -< 
ip*{X) is the same as Lp*{w~^ii) -< Lp*{w^^X). This imphes that w^^fi <" w^^X, and hence 
fi < X. The last statement follows directly from the second. □ 

Theorem 4 Let (/^i , : g — )■ g &e two embeddings of q into g. Let Xi,yi, hi for 1 < i < r 
form a canonical set of generators of g. Write x},yl,h} and xf.yf.hf for their images 
under tpi, respectively. Assume that h},hf G f). Then ipi and (^2 o'^e linearly equivalent 
if and only if there is a w E W with w{h\) = hf for 1 < i < r. 

Proof. First suppose that w E W exists. Let p : g — )■ g[([/) be a representation of 
g. Let /X G be a weight of p, i.e., there are nonzero u E U with p{h)u = p{h)u for 
all /i G f). Set Pi = pLpi for i = 1,2. Then pi is a representation of g. Observe that 
p{hj) = K{p,h}) = K{wp,w{hj)) = K{wp,w{h})) = {wp){w{h})). But also wp is a weight 
of p, with the same multiplicity. Hence it follows that pi and p2 have the same weights 
with the same multiplicities. Hence ipi, <^2 are linearly equivalent. 

Now assume that </?i, v?2 are linearly equivalent. By Lemma |3] there are 0^,02 G A^gC)) 
such that = aiipi have the properties stated in Lemma [3] for i/j. 

Let /9 : g — )■ Qi{U) be an irreducible representation with highest weight A. Set pi = poipi. 
Then from Lemma |3] it follows that '0*(A) is the largest weight of pi in the ordering -<. 
As the are linearly equivalent it follows that (A) = "02 (-^)- This is the same as 
X^t/ji^h)) = X{ip2{h)) for all h E i). Now since f)^ is spanned by dominant weights, this 
equality follows for all A G f)^- Hence i/Ji{h) = ■02(^) for all h E I). In particular, this is 
true for the hi. So ai{(pi{hi)) = cr2(v?2(^j))- Now let Wi E W he such that (Tj|(, = Wi. Then 
we get the statement of the theorem with w = W2^wi. □ 

Corollary 5 Let Qi, Q2 be two semisimple subalgebras of g, both isomorphic to g. Let 
h\, . . . , hi, hi, ... , hi, be the h-parts of canonical sets of generators o/gi andQ2 respectively. 
Assume that h^ E i) for all i, k. Then Qi, Q2 o'^e linearly equivalent if and only if there is 
a w E W with 

{w{hl) \ I < i < r} = \ 1 < i < r}. 

Proof. The "if" -part follows immediately from Theorem |H For the "only if " -part suppose 
that gi, g2 are linearly equivalent. Let p : q sK^) be a faithful representation. Then 
there is a G GL{U) with ap{Qi)a~^ = p{q2)- Set hf = p~^{ap{h])a~^) E Q2- Let G2 C G 
be the inner automorphism group of g2. Let [)|, i)l denote the subspaces of g2 spanned 
respectively by the h^ and the hf. These are Cartan subalgebras of g2 so there is a a G G2 
with cr{i)l) = f)2. Set hf = cr{h\). Then also the hj form the /i-part of a canonical set 
of generators of g2, lying in the same Cartan subalgebra of g2 as the hf. Let 1^2 denote 
the Weyl group of g2 with respect to \)2- Since different sets of simple roots of g2 are 
conjugate under W2, there is a m G W2 such that {u{h1)} = {h^}. Let r E G2 he such 
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that r restricted to f)2 is u. Let hi, . . . ,hr form the /i-part of a canonical generating set 
of Q. Let : — )■ 01 be the isomorphism sending /ij to hj. Set (y92 = ratpi. Then is 
hnearly equivalent to ^pi. Moreover, {ip2{hi)} = {hf}. Now we get the required w E W 
from Theorem HI □ 

So we can decide linear equivalence if we can decide whether two sets of elements, 
{h\, . . . , h}.} and {hf, . . . , h'^} of f)ffi are conjugate under W. Since W preserves the Killing 
form we assume that the ordering is such that K{hl, h^) = K,{hf, h'j) for 1 < i, j < r. 

We fix a set of positive roots and corresponding set of simple roots {ai, . . . ,«/}. 
Then the reflections Sq,. generate W. 

From Section [ZD we recall that C C [)m is the set of all h with a{h) > for all a > 0. 
We note that for a given h G [)^ it is straightforward to find its unique VF-conjugate lying 
in C. Indeed, initially we set Hq = h. Let i > and suppose that hi is found. If /ij G C 
then we are done. Otherwise there is aj with aj{hi) < 0. Then set /ij+i = So,.{hi). Note 
that for i < j we have hj — hi = ^'^^^ flfcOfc with G M non-negative, and at least one 
coefficient ak is positive. Hence all hi are different, and as W is finite the sequence of the 
hi must land in C. From this we also immediately get a w & W with w{h) G C. 

Next we have a method for deciding whether there is a. w & W with w{h\) = h^. We 
first compute tf 1,102 G W with Wi{h\) G C. If those are not equal, then the required w 
does not exist. Otherwise set u = W2^wi; then u{h\) = h\. Now the set of all t; G 
sending h\ to h\ is exactly '^tahw{h\)u, where StabvK(^i) denotes the stabiliser of h\ in 
W. 

Set h = W2{hl). Let / be the set of all i with ai{h) = 0. It is known (cf. [9J, Theorem 
1.12) that Stabw{h) is generated by the Sq- with i E I. Now Stabvi/(^i) = W2^Stabw{h)w2- 
This implies that Stabwihf) is generated by the refiections s^-i^^.-j, where z G /. The roots 
w^^(aj) for 2 G / form a simple system of a root subsystem of $, of which Stahwihl) is 
the Weyl group. 

Now set /if = ^(/li) for 1 < i < r. We decide if there is i; G Stahw{hD such that 
v{hi) = hf for 2 < i < r. We can do this as the sequence is shorter. If such a v exists, also 
the required w (which is vu) exists. In the other case it does not. 

Finally, in order to decide whether there is a to G with {w{h})} = {hf} we loop over 
all permutations tt of {1, . . . , r} with ^{hl^^-^, h\f^--^) = K.{h}, /ij) for 1 < i, j < r. For each 
such 71 we decide whether there is a w E W with w{hl^^-^) = hf. Once we find one we stop. 
Remark. This procedure works well in partice if the number of permutations as above is 
small. This very often is the case. The main exception being the case where {h\, . . . , h]) 
is the /i-part of a canonical generating set of a Lie algebra of type kAi. In situations like 
that the algorithm has to work a lot harder, as up to k\ permutations have to be tried. 
Fortunately, for the simple Lie algebras of ranks < 8 there are not many subalgebras of 
such a type with large k. 

Remark. If a class of linearly equivalent subalgebras splits into more than one class of 
equivalent subalgebras, then each of the latter classes has a representative having a canoni- 
cal generating set with /i-part that is the same for each of them. Only the other generators 
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Xi, Hi differ. 



4 Constructing a subalgebra 

In this section we describe algorithms for constructing a canonical generating set of a 
semisimple subalgebra of q, given its Cartan matrix and /i-part. 
Let [) be a subalgebra of [). For /i G [)* we set 

= {x e Q \ [h,x] = fi{h)x for all h Ei)}. 

Then g is the direct sum of the various 

Lemma 6 Let E i)* be such that g{fi) ^ and such that there is a h E\) with fi{h) = 2. 
Set 

0^ = {ueg{ti) I [0(O),^i] = g(^)}, 

= {e E g(/i) I 3/ G 0(— /i) wi/i {h, e, /) zs an sl2-triple}. 

Let Go be the connected subgroup of G with Lie algebra 0(0). Then Gq has a dense orbit 
in 0(/x), which is equal to O^. If is nonempty then E^ = O^. 

Proof. By standard arguments it is proved that k, is non-degenerate on 0(— /i) ©0(/i), and 
on 0(0). Hence 

fcez 

is a reductive Z-graded Lie algebra. In [17] it is shown that 0(^) has a dense Go-orbit. It 
is clear that a u G 0(^) lies in this dense orbit if and only if it lies in (9^. 

Suppose that i?^ is not empty, and let e G i?^. Then from s[2-representation theory it 
follows that ade : 0(0) 0(/i) is surjective. In other words, [0(0), e] = 0(/i). Hence the 
Go-orbit of e is dense in 0(/i). So this last orbit coincides with O^. But then also E = O^. 

□ 

Let C be the Cartan matrix of the root system of a semisimple Lie algebra. Let 
hi, . . . ,hr G f). We want to find Xi,yi G satisfying the relations ([I]), or decide that no 
such elements exist. We assume that hi eH, as otherwise the required Xi,yi certainly do 
not exist. The space spanned by hi, hr will be^enoted f). 

First of all, let /ij G f)* be defined by ^iihj) = C{i,j). We compute bases of 0(/Xj) and 
0(— /ij), and of 0(0), which is the centralizer of f). The Xi, jji, if they exist, lie in 0(//i), 
0(— /ij) respectively. 

In the second step we find xi G 0(/ii), jji G 0(— /ii) such that {hi,xi, yi) is an s[2-triple. 
For this we use Lemma [61 After trying a few random elements we find an xi G 0(/ii) with 
[0(0), = 0(yUi), i.e., such that xi lies in (9^^. By solving a set of linear equations we 
either find yi G 0(— /ii) such that {hi, Xi,yi) is an sl2-triple, or we decide that that no such 
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yi exists. In the latter case there is no sl2-triple {hi,xi,yi) with xi G yi G g(— /ii)- 

Indeed, in that case the set E^^ (notation as in Lemma [6]) is empty. So in the latter case 
we stop with the conclusion that the Xi, yi do not exist. In the former case we continue. 

In this second step we choose a random element xi. We do stress that for the existence 
of the subsequent elements for z > 1 it does not matter which xi is chosen, as long 

as [g(0),xi] = Indeed: all elements with that property are conjugate under G(0) as 

they lie in the same dense orbit. 

Now we continue to find the remaining Xj, yi. For this we use two methods, which we 
call the linear method and the polynomial method. 

For the linear method we suppose that Xi,yi, 1 < z < s, for a certain s with 1 < s < r, 
have been found, satisfying ([1]). We also assume that all different such sets are G-conjugate. 
In other words, if x^, y[ for 1 < z < s also satisfy ([1]), then there exists cr G G with (j{x[) = Xj, 
(t(?/^) = jji, a{hi) = hi. Note that by the above construction this is certainly true for s = 1. 

Set 

= {uE sifJ-s+i) I [u, i/i] = for 1 < i < s}, 

g\-^s+i) = {u e Q{-fis+i) I [u, Xi] = for 1 < i < s}. 

Then Xg+i G g'{iJ,s+i), Vs+i £ q'O^s+i)- Let also q'{0) be the intersection of 0(0) and the 
centralizer of all Xi, 1 < i < s. By s[2-representation theory it follows that q'{0) centralises 
also all jji, I < i < s. Hence 0'(O) acts on ^'(/is+i). Let Gq be the connected subgroup of 
G with Lie algebra q'{0). There are now two cases that can occur. 

In the first case, after trying a few random elements, we find a x^+i G g'^j-s+i) with 
[^'(0), Xs+i] = Q'{fis+i)- This means that Gq has a dense orbit in 0'(/is+i)- By solving 
a set of linear equations we either find jjs+i G Q'{fis+i) such that {hs^i,Xs+i,ys+i) is an 
sl2-triple, or that no such ijs+i exists. In the former case we say that the linear method 
has successfully found Xg+i, jjs+i- Note that this also implies that all sets of Xj, jji for 
1 < i < s + 1 are G-conjugate. In the latter case we say that the linear method has broken 
down at step s + 1. 

The second case occurs when, after trying a few random elements, we do not find an 
Xg^\ as above. In this case we also say that the linear method has broken down at step 
s + 1. 

After having found Xi, yi we repeat the linear method. If it does not break down then 
in the end we find a complete set of Xi, yi. If it breaks down at step s + then we use the 
polynomial method. 

So for the polynomial method we also assume that Xij^j, for 1 < z < s, have been 
found, satisfying the relations ([T]). For s + l<A;<rwe compute bases of the spaces 

Q'{^^k) = {m G s(/ifc) I [u, i/i] = for 1 < i < s}, 

0'(-/«fc) = {ue Q{nk) I [u, Xi] = for 1 < z < s}. 

We express the Xj, fji for s + l<z<ras linear combinations of the bases of, respectively, 
g'(//fc) and 0'(— /Lt^), with indeterminates as coefficients. Then the Xj, jji satisfy ([T]) if and 
only if certain polynomial equations in the coefficients are satisfied. We compute the 
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polynomial equations, and by Grobner basis techniques (see Section [215]) . we either solve 
them, or decide that no solution exists. 

Remark. Note that the linear method is heuristic in nature. However, it is automatic: 
if it succeeds then no further intervention is necessary to construct the subalgebra. We 
note also that there are situations where the linear method must break down as there are 
subalgebras that are only defined over an algebraic extension of Q. In this case using the 
polynomial method is necessary. However, this last method is not entirely automatic (cf. 
Section [231) . 

Remark. In the next section we give methods to construct a suitable set of candidates 
{hi, . . . , hr) for the /i-parts of canonical generating sets of semisimple subalgebras. This 
construction is such that {hi, . . . , /ir—i) will be the h-part of a canonical generating set of 
a subalgebra of rank r — 1. However, the Xi, yi for l<z<r — Ido not necessarily lie in 
the bigger subalgebra, as in the two cases (the algebra of rank r — 1 and of rank r) the 
spaces 0(;Uj) are quite different. 



5 Finding candidates 

In this section we deal with the problem of finding a suitable set of candidates for the h- 
parts of canonical generating sets of semisimple subalgebras of g, with given Cartan matrix 
C. For this we first consider a problem involving characters, whose solution will help us in 
making the set of candidates smaller. 

Let be a semisimple Lie algebra with canonical generators Xi, jji, hi, 1 < i < r 
satisfying ([T]). Let be a finite-dimensional g- module. Then V is spanned by common 
eigenvectors of the hi. Moreover, the eigenvalues of the hi are integers. For an e = 
(ei, . . . , e^) G Z'' we set 

Ve = {v & V \ hi ■ V = CiV for 1 < z < r}. 
Let be indeterminates and write I'' . Then the polynomial 

^(dim\4)x" 

eeZ'- 

is called the character of the g-module V . 

Now fix i with 1 < i < r. For m G Z we set = {v (zV \hi-v = mv}, and define the 
polynomial 

/.(x,) = 5^(dimVjxr. (2) 

We call the polynomial /i + ■ ■ ■ + /r the character-puzzle of V . It is clear that from 
the character of V we can compute its character-puzzle. More generally we say that a 
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polynomial of the form fi{xi) + - ■ ■+fr{ character-puzzle. It is clear that a character- 

puzzle does not necessarily correspond to a character. If it does we say that it is solvable. 
Here we consider the following problem: given a character-puzzle / = fi{xi) + ■ ■ ■ + fr{xr) 
decide whether it is solvable. 

For this we proceed as follows. First we note that ^ is a direct sum of simple modules, 
determined by a highest weight, which is an e = (ei,...,e.r) with Cj > 0. From the 
character-puzzle we retrieve all non-negative eigenvalues of the hi. This gives a finite 
number of possibilities for the highest weight of a simple constituent of V. For each possible 
highest weight we compute the character of the corresponding highest weight module (cf. 
[7J), and from that its character-puzzle g. Then we subtract, h = f — g. Then recursively 
we establish whether h is solvable. 

If at least one h that we so obtain is solvable then / itself is solvable. Otherwise it is 
not. 

Let C be the Cartan matrix of (the root system of) a semisimple Lie algebra q of rank 
r. In this section we describe how we find a set H of r-tuples {hi, ... , hr) G W such that 
every semisimple subalgebra of g isomorphic to q is linearly equivalent to a subalgebra 
with canonical set of generators Xj, y^, hi with {hi, . . . ,hr) G H. We also want the set 
to be "small" (whatever that means). So, although the set would be a solution to 
the problem, it is far too big. (For example, if q is of type Es then it has 2611951200'' 
elements.) 

A first reduction is given by CorollaryO if there are two r-tuples {hi, ... , hr), {h'l, . . . , h'^.) 
such that there is a w & W with w{hi} = {h'i}, then we can discard one of them. 

Secondly, let Cq be the (r — 1) x (r — l)-matrix in the top left corner of C. Then we 
may assume that we know a set Hq of (r — l)-tuples {hi, ... , hj—i) G such that every 
semisimple subalgebra of g with Cartan matrix Cq is linearly equivalent to exactly one 
subalgebra with canonical set of generators Xj, jji, hi, 1 < i < r — 1, with {hi, ... , hr-i) G 

Therefore we only put r-tuples {hi, ... , hr) into the set H that have {hi, . . . , hr-i) G Hq. 
Note that for r = 2 we know the set Hq from the classification of the nilpotent orbits in q. 

So let {hi, ... , hr-i) G Hq. We want to extend this (r — l)-tuple with an hr. If we just 
take any hr then the set H gets too big. So we perform further reductions. For this 
we distinguish two cases. ^ 

In the first case, in the Dynkin diagram of C, the node labeled r is not isolated. So it 
is connected with 1, 2, or 3 bonds to a simple component Fq of the Dynkin diagram of Cq. 
Let F be the simple component of the Dynkin diagram of C, containing Fq. Let ii, . . . ,is 
be the labels of F, where is = r. Let g be a simple Lie algebra with Dynkin diagram F, set 
of canonical generators Xi, iji, hi, for 1 <i < s, and Killing form k. As seen in Section 1^751 
the matrix {K{hi^, hi^)) is a scalar multiple of the matrix B = {k{hi, hj)). Furthermore, we 
know the scalar factor r] from comparing nijii^, hi^) and k{hi, hi). In particular we know 
what K{hr, hr) has to be; denote this value by 6. 

Now let hi, . . . ,ht be representatives of the ly-orbits in 7/ (see Section [2^ . Note that 
k{u,u) = K{hi, hi) for all u in the W^-orbit of hi. So we enumerate the orbits of those hi 
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such that K{hi, hi) = 9. A hr in such an orbit is selected if the matrix {K{hi^, hi^)) is equal 
to 9 times B, and K{hr, hi) = for i not in {ii, . . . ,is} (cf. Lemma [1]). 

If the number of bonds is 1 then we can reduce the work further. Suppose that the 
node labeled = r is connected to the node with label is-i in T. Let . . . , /3s be the 
simple roots of q. Then f3s-i and are conjugate under the Weyl group W of q. Also, 
/3/ = hi. So from what is said in Section 12.11 it follows that hg-i and hs are conjugate 
under W. Hence they are conjugate under G, the inner automorphism group of g. Now 
an embedding g g induces an embedding G ^ G. It follows that hi^_^ and hj, must 
be conjugate under G, which implies that they are conjugate under W . The conclusion is 
that we can limit our search for suitable elements hr to the W^-orbit of hi^_^. 

In the second case, in the Dynkin diagram of C, the node labeled r is isolated. In other 
words, a subalgebra isomorphic to q is the direct sum of a subalgebra Qq, with Cartan 
matrix Cq, and a subalgebra isomorphic to s[2. Then by Lemma [H we can restrict to 
adding the hr with hr) = for 1 < i < r — 1. Also in this case we run through "H by 
enumerating the VF-orbits of the hi. 

In both cases we can still encounter VF-orbits that are too large to enumerate. For 
example, in order to construct the subalgebras of type A2, or of type 2^41, with the above 
procedure, one would have to run through all orbits; when g is of type Eg this amounts 
to examining 2611951200 elements. In order to reduce the work needed we use character- 
puzzles. Let V be the smallest nonzero g- module. For {hi, . . . , hr) G Ti^ we compute the 
corresponding character-puzzle, as in ([2]), where we view V a.s a g-module. We note that 
all hr in the VF-orbit of hi lead to the same character-puzzle. So we decide if the character- 
puzzle of corresponding to {hi, . . . , hr-i, hi) is solvable (i.e., corresponds to a character of 
g). Only if it is, we enumerate the orbit of hi. 

Remark. The procedure using character-puzzles eliminates the largest orbits. For exam- 
ple, for g of type Es, there are 11 orbits (out of a possible 69) that need to be enumerated 
for constructing the subalgebras of type A2; they have sizes 240, 2160, 6720, 17280, 30240, 
60480, 69120, 181440, 241920, 483840, 1814400. We also note that Eg is a difficult case 
in two respects: it has by far the largest Weyl group, and the largest minimal faithful 
representation of all simple Lie algebras of ranks < 8. The fact that the minimal faithful 
module has dimension 248 makes solving the character puzzles rather hard. However, it is 
still worth the wile, as the orbits that are excluded this way are so big. From the sizes of 
the orbits that still need to be enumerated we also see the need for an algorithm, as the 
one of Snow ([IS]), that does so using little memory. 



6 Subalgebras of semisimple Lie algebras 

Let g = gi © g2 be the direct sum of two semisimple ideals. Let f) = f)i © f)2 be the 
corresponding decomposition of the Cartan subalgebra. Then the Weyl group of g is 
a direct product Wi x W2, where Wi (respectively W2) acts trivially in f)2 (respectively 
P)i). Let Li be the set of representatives of the hnear equivalence classes of semisimple 
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subalgebras of 0j. We assume that each element of Li has a canonical set of generators 
with h-pait lying in i). 

Let a© bi, 62 © c be elements of Li, L2 respectively, where bi, b2 are isomorphic. Let 
h},x},yl for 1 < i < s, h},x},y} for s + 1 < i < s + r, hf^xl^yf for 1 < z < r, h^,x^,yf 
forr + l<z<r + mbea canonical generating sets of respectively a, bi, b2, c. We 
assume that the canonical generators of bi, b2 are "in the same order"; that is, mapping 
— )■ x], —J- yf, h\j^^ hf, ior 1 < i < r defines an isomorphism bi — )■ b2. Let ir be 
a permutation of {1, . . . , r} preserving the Cartan matrix of b2, or, equivalently, such that 
^2(^^(4), ^^(j)) = ^2(/if, h'j) for I < i, j < r, where K2 denotes the Killing form of 02- Then 
also mapping xl_^i x^^.^, yl^^^ j/^^.^, /i^.^^ /i^^.-, defines an isomorphism bi b2. Let 
now Q be the subalgebra of g with canonical generating set 

{h],x],yl \ l<i<s}U {hl^^ + x^^^ + 1/]+^ + | 1 < i < r}U 

{h'^,x'^,yf \ r + 1 < i < r + m}. 

Let L denote the set of subalgebras of q that can be constructed this way. The next 
proposition is similar to [5j, Theorem 15.1. 

Proposition 7 Every semisimple subalgebra of g is linearly equivalent to an algebra in L. 

Proof. Let g be a semisimple subalgebra of g. We may assume that it has a canonical set 
of generators with /i-part lying in f). 

Let Pi : g ^ Qi denote the projection homomorphism. Then kerpi fl kerp2 = 0. Let a, 
c be the sum of the ideals of g that lie respectively in kerp2 and in kerpi. Let b be the 
sum of the remaining ideals. Then g = a © b © c, with a C gi, and c C g2. Let h\, . . . ,hl, 
h'^^i, . . . , /i^+m denote the /i-parts of canonical generating sets of a and c respectively. Then 
h} G /i^ e fj2. 

Let hi,Xi,yi, 1 < i < r be a canonical generating set of b. Let H denote the set 
containing the hj, 1 < i < s, h"^, r + 1 < i < r + m and hi, 1 < i < r. We must show that 
there is a. w E W such that w{H) is the /i-part of a subalgebra in L. 

Note that pi and p2 are injective on b. Write hl^i = pi{hi), /if = P2{hi), for 1 < i < r. 
Those elements form the h-part of a canonical generating set of a semisimple subalgebra 
bi, respectively b2, of gi and g2. Moreover, the b, are isomorphic to b. In particular, 
h\, . . . , hl_^_J. form the /i-part of a canonical generating set of the subalgebra a © bi of gi. 
Therefore, after possibly reordering the elements of H, there is a G Wi such that the 
Wi{hj) form the /i-part of a canonical generating set of an element of Li. Note that this 
fixes the ordering of the hi G H. We can still reorder the hf G H, where r + l<i<r + m. 
So there is a 1^2 G W2 such that W2{h\), . . . , tf2(/if), '«^2(^f+i)> • • • )'^2(^f+m) form the h- 
part of an element of L2, up to, possibly, a permutation vr of the first r elements. This 
permutation has to leave the Cartan matrix of b2 invariant. □ 

Proposition [7] gives an immediate procedure for finding a set L containing representa- 
tives of all linear equivalence classes of semisimple subalgebras of g. However, it can still 
happen that different members of L are linearly equivalent. For weeding out linear equiv- 
alent pairs we use the algorithm outlined in Section |3l We also note that L is the disjoint 
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union of two subsets L', L" . Here V contains the subalgebras that are the direct sum of 
an algebra in Li and an algebra in L2. And L" has the algebras constructed as above with 
^'1,^2 7^ 0. Among the algebras in L' there are no linear equivalences. Furthermore, an 
algebra in V is never linearly equivalent to an algebra in L" . 

7 Deciding inclusion 

For two semisimple subalgebras fli, 02 C we write gi — 02 if 0i is linearly equivalent to a 
subalgebra of 02- (Here linear equivalence is defined with respecto to 0.) Given 0i, 02, with 
canonical generating sets with /i-parts in f), we decide whether 0i — )■ 02 in the following way: 
First we let L be the set of representatives of the classes of linear equivalent subalgebras of 
02- We get this from the classification of those subalgebras of 02- All are given by canonical 
generating sets having /i-parts in f). Then we decide whether 0i is linearly equivalent to 
an element of L, using the algorithm from Section Ei 

Now let 01, . . . , 0s be a chain of subalgebras. This means that Qi — t- 0^+1 for 1 < z < s. 
Then we can compute a realization of the chain; that is if necessary we replace the Qi by 
linear conjugates such that Qi C 0i+i. For this we start "at the top", and suppose that 
0j — 7- ■ • ■ — 7- 0s has been realised. We compute the subalgebras of Qi isomorphic to 0i-i, 
up to linear equivalence. We find a subalgebra s that is linearly equivalent to 0i-i, and 
replace 0j_i by s. 
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